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We explore the possibility that weakly interacting dark matter can form bound states - WIMPo- 
nium. Such states are expected in a wide class of models of particle dark matter, including some 
limits of the Minimal Supersymmetric Standard Model. We examine the conditions under which we 
expect bound states to occur, and use analogues of NRQCD applied to heavy quarkonia to provide 
estimates for their properties, including couphngs to the Standard Model. We further find that it 
may be possible to produce WIMPonium at the LHC, and explore the properties of the WIMP that 
can be inferred from measurements of the WIMPonium states. 

PACS numbers: 14.80.-j 

I. INTRODUCTION 

The identity of the missing mass of the Universe is one of the most pressing questions at the interface of 
cosmology and particle physics. The existence of dark matter is one of the few noncontroversial indications 
for physics beyond the Standard Model (SM) , with strong preference for a weakly interacting massive particle 
(WIMP) as the best explanation for all of the evidence. The hunt is on to discover WIMPs through indirect 
detection of its annihilation products, direct detection of its scattering with nuclei, and production at high 
energy accelerators. In particular, the seeming coincidence between WIMP masses at the TeV scale and 
the scale of electroweak symmetry breaking would seem to argue for a connection between the two. That 
connection, in turn, would suggest that WIMPs should not interact only via gravity, but could have large 
couplings with Standard Model fields. The realization that the thermal prediction for the cosmological relic 
density of WIMPs would roughly explain the observations lends further credence to the scenario. 

While a successful WIMP should be weakly interacting with the SM, constraints on the interactions 
between the WIMPs themselves (or among particles of the entire "dark" sector) are quite weak. The relic 
density suggests a size for the interaction between WIMP and ordinary SM particles, but says little about 
the interactions among the WIMPs themselves. Quite generally, there may be large forces between WIMPs 
mediated by light particles, and their exchange can induce a Yukawa potential between WIMPs, which may 
be attractive. In some cases, there will be bound states of WIMPs as a result - "WIMPonium" . When 
WIMPs are literally weakly interacting (that is, charged under the SU{2)l x U{1)y of the Standard Model), 
they will often couple to the Z boson, and perhaps also to the Higgs. WIMPonium is potentially a generic 
feature of models of dark matter. 

In this article we examine the circumstances under which we expect WIMPonium to have interesting 
consequences, and explore some of the phenomena which are expected to result. Early studies have focused 
on the specific scenario of WIMPs charged under the electroweak SU{2) x U(l) and examined the effects on 
the prediction for indirect detection signals [T] or the relic density [5]. It has also been recently suggested 
that WIMPonium states, either realized explicitly or manifest as nearby virtual states in the form of a large 
Sommerfeld effect, can reconcile observations by PAMELA and ATIC [T with a thermal relic WIMP [HEHB] 
(however, see [7j). We take a somewhat more generic viewpoint, and allow for general interactions between 
the WIMPs without recourse to a very specific model of dark matter. We also (for simphcity) restrict 
ourselves to the case when the WIMP's SU{2) partners (if any) are not highly degenerate in mass with the 
WIMP, and consider effective theories containing only the WIMP itself. 

This article is organized as follows. In Section|ll] we review the properties of the bound states of the Yukawa 
potential, and discuss the conditions under which we expect WIMPonium states to occur in the spectrum. 
In Section |III| we make rough estimates for WIMPonium production and decays, using the language of 
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FIG. 1: The minimum value of for which there is an Is WIMPonium state (sohd hues), 2s states (dashed lines), 
and 2p states (dotted lines) as a function of the WIMP mass M and for Yukawa screening mass (bottom curve to 
top curve) m = 1, 10, and 100 GeV. 



NRQCD ^ as previously applied to systems of heavy quarkonia [9]. Section IV considers the possibility 
that WIMPonium states could be produced at future colliders and what such states could reveal about the 
underlying theory of dark matter, and Section [V| explores their impact on astrophysical observations. We 
conclude in Section IVll 



II. WIMPONIUM CHARACTERISTICS 



The Yukawa potential between two WIMPs is characterized by its strength and range m, 

V{r) = -^e-™"" , (1) 

where our convention is such that > corresponds to an attractive potential, and < is repulsive. 
The nature (attractive or repulsive) of the potential acting on a pair of WIMPs depends on the spin of 
the WIMP and that of the exchanged boson (which we will refer to as the "mediator" particle, generically 
denoted by t^). As is well known, a scalar exchange leads to a universally attractive potential, whereas a 
vector exchange will attract particles to antiparticles, but repel pairs of particles or pairs of antiparticles. 

For a weakly coupled system, the parameter can be computed from the non-relativistic limit of a single 
particle exchange (see, e.g. |TD] for a pedagogical review). For complex scalar or Dirac fermion WIMPs 
exchanging a vector, = /^tt where g is the gauge coupling of the exchanged vector, {g = \/ g\+ g2 for 
a Z-boson exchange). Similarly, a fermion WIMP exchanging a scalar mediator will have given by the 
strength of its Yukawa interaction with the scalar, = |yp/47r. A scalar WIMP interacting with a scalar 
mediator via a three-point interaction has a coupling A with dimensions of mass in natural units, and the 
strength of the induced Yukawa potential is given by = |Ap/(47rM2) where M is the mass of the WIMP. 

Unlike the Coulomb potential, a non-relativistic system governed by the Yukawa potential has a finite 
spectrum of bound states, whose energy levels depend on the orbital angular momentum I. Detailed numerical 
solutions to the Schrodinger equation for the Yukawa potential are presented in [TT], and we describe some 
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general features of the solutions here. The characteristic size of the system can be expressed in terms of the 
corresponding Bohr radius for two WIMPs of mass M, 

a- = ^ . (2) 
The system can be characterized by how "Coulomb-like" it is through the quantity, 

D=^ = ^. (3) 

am 2m 

In the limit D ^ \ the lowest lying states have wave functions small compared to the radius at which the 
exponential suppression of the Yukawa potential is noticeable, and their properties are very similar to those 
of the Coulomb potential. The quantity D also controls the spectrum of bound states. In order to have a 
bound state at all (the Is state), the potential must satisfy D > 0.8. The 2s state requires D > 3.2, and the 
2p state will be present provided D > 4.5. In figure [l] we present the minimum for which Is, 2s, and 2p 
bound states exist as a function of the WIMP mass M for several values of mediator mass m. 
The mass of the WIMPonium state is — 2M — Eb where Eg is the binding energy, 

alM 

EeinJ) = -^e{n,l) . (4) 

The function e{n,l) asymptotes to in the limit D ^ 1. For the Is state, e = {0.02,0.1,0.8} for 

D = {1, 1.4, 10}. As expected in a weakly bound system, the binding energy is typically quite small compared 
to the mass of the WIMPs themselves; for a WIMP of mass 1 TeV with ~ 10"^ and D ^ 1, the binding 
energy of the Is state will be Eb ~ 25 MeV. 

While limits on are very weak or non-existent, we do expect WIMPs to be reasonably weakly coupled, 
and we can generically expect that there will be a limited number of WIMPonium states. The lowest states 
are the Is, the 2s, and the 2p states. Since the Is and 2s states possesses no orbital angular momentum, the 
spin of the WIMPonium will be determined entirely by the spins of the WIMPs themselves. As symmetric 
states, they further require that the spins of the WIMPs enforce either Bose-Einstein or Fermi-Dirac statistics. 
Consequently, scalar WIMPs will result in scalar s-wave WIMPonium states. Majorana fermions are required 
by their statistics to be in a scalar spin state. Dirac fermions could potentially form either particle-particle 
(antiparticle-antiparticle) bound states (which would be required to be scalar) or particle-antiparticle bound 
states, which could be scalars or vectors. A 2p state would allow a complex scalar to form a particle- 
antiparticle bound state which is a vector. We will often use the usual spectroscopic notation "^^^^Lj to 
label the angular momentum configurations of the lowest lying bound states. 



III. WIMPONIUM EFFECTIVE THEORY 



We compute the rates for WIMPonium production and decay by using non-rclativistic quantum field 
theory, borrowing well known results from the study of heavy quarkonium [9, . We make use of the fact that 
for weak coupling, ^ 1, the velocity of a WIMP inside a WIMPonium state is of order v ^ a^, with 
good separation between the rest energy M, the typical WIMP momentum, Mv ^ Ma^, and its kinetic 
energy, Mv'^ ~ Ma^. WIMPonium properties may thus be thought of as a perturbative expansion in both 
couplings and velocity. 



A. WIMP Interactions with the Standard Model 



Before considering the effective field theory for energies much less than M, we consider the relativistic 
effective theory at scales of order M which describes the WIMP interactions with SM fields. These inter- 
actions typically appear as higher dimensional operators involving a pair of WIMPs (to insure the WIMP 
is stable with respect to decays into SM states) and a local product of SM fields. Such contact interactions 
are typically generated from a UV complete theory by integrating out heavier degrees of freedom. For ex- 
ample, in the MSSM a neutralino-neutralino-e-e interaction arises from integrating out the selectrons. As 
long as the masses of these additional states are somewhat heavier than that of the WIMPs themselves, this 
effective field theory should be a reasonably accurate description of the SM interacting with pairs WIMPs. 



3 



It is relatively straight-forward to generalize our results here to the case in which there are additional dark 
sector states close in mass to the WIMP (as will be inevitable, e.g. in a scenario employing co-annihilation 
to arrange for the correct relic density) by explicitly including such states in the effective theory. 

A pair of complex scalar WIMPs can interact with Standard Model fields through operators (up to 
dimension six), 

f K f ) 

+ ^\X?\H\'' + (x*^,x) {H^D^H) + -^W^W.^W^'^ + -^Ixl'B^^B^'^ + H.c. (5) 

where the parameters and A represent combinations of masses and couplings in the UV theory, and H is 
the SM Higgs doublet, including the would-be Goldstone bosons. Some of these interactions were previously 
considered (in a different context) in (12j. We have assumed that there are no new sources of chiral flavor 
symmetry breaking, and thus the operators containing fj^fa are flavor conserving and accompanied by the 
fermion Yukawa interaction y/. We have assumed the complex scalar WIMP has a conserved global U{1)^ 
symmetry (which could be invoked to motivate the fact that the WIMP stable). If such a symmetry is 
absent, one can also write terms (with independent coefficients) with |xP ^ X^- A real scalar WIMP has 
only such terms, and in addition, the terms containing x*9^X m^^Y be expressed in terms of other terms 
already present by integrating by parts and applying the equations of motion [T3]. 

A Dirac WIMP has potentially both scalar and vector interactions with SM flelds (up to operators of 
dimension seven), 

^ = E I [X^'^X] UliM + ^ [xrx] Ub^iM + ^ [XX] HJM + ^ m \H\' 

f { fh /h / J ^ 

+ ix^x] {H^D^H) + [xx] \H\' + [x^d'^x] B,. + [xx] W^.W^^ 

+ ^[XX]B^.B'^'^ + H.c. (6) 

^^B 

There are also pseudoscalar (x75X) or pseudovector (xi^lsx) versions of each operator. A Majorana WIMP 
expressed as a four component spinor looks much the same, though the X7^x structure vanishes when the 
WIMP is Majorana. 

B. Non-relativistic Effective Theory 

To construct the non-relativistic effective theory, we express the WIMP in terms of a fleld which annihilates 
particles and one which creates antiparticlcs, with the large time oscillatory behavior due to the mass factored 
out, 

X - ^[e-^"e + e^*"'7] (7) 

(8) 



X 




2M 



piMt _ ■ -iMtS-V c 



where ^ annihilates a WIMP and 77 creates an anti-WIMP, and the two lines refer to a scalar or fermion 
WIMP, respectively. In the fermion WIMP case the non-relativistic flelds are two component Pauli spinors 
and (7 are the Pauli matrices. For a real scalar WIMP -q ^*/\/2 (and ^ — > C/\/2) and for a Majorana 
fermion, r/ C^/v^ = — zcr'^^*/\/2 (and ^ — > ^/\/2). These decompositions lead to kinetic terms in the NR 
limit in which the fast oscillations have been explicitly removed. 



2M J 



V (9) 



with some mild abuse of notation since for a scalar WIMP these are scalar flelds, and for a fermion they are 
Pauli spinor flelds. Higher time derivatives are explicitly removed by use of the equations of motion, and 
result in higher order (in 1 /M) relativistic corrections to propagation ^ . 
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WIMPonium production and decay are both captured by the imaginary coefficients of local four- WIMP 
operators which arise in the UV theory as "box" diagrams in which two WIMPs annihilate into SM fields 
which then recombine into two WIMPs. In terms of the couplings defined in equations (jsj and ^ they 
generally arise as two applications of a given term with the Standard Model fields contracted into a loop. 

For a scalar WIMP, the relevant operators for WIMPonium decay take the form, 

Jf2 U{So) + 0{ Fi), (lUj 

i 

where we have included the intermediate states //, H*H, WW, and BB. The four-WIMP operators 
encoding WIMPonium decays are, 

o+CSo) = {CC)\o){o\m, (11) 

0-CS„) = {VV)\0){0\{V*V*), (12) 
O'CSo) = (^r)|0)(0|(,7*0, (13) 
O(iPi) = (r7*Ve)|0)(0|(rV77), (14) 

corresponding to s-wave particle-particle, antiparticle-antiparticle, and particle-antiparticle states and a p 
wave particle-antiparticle state, respectively. The Wilson coefficients Ci may be computed in terms of the 
interactions of Eq. ([5| and will have imaginary parts (through the optical theorem) when a pair of WIMPs is 
able to non-relativistically annihilate into the given channel. There will also often be annihilations into the 
light particle responsible for binding the WIMPonium system which may make up an appreciable branching 
fraction. 

WIMPonium production is described by similar operators. For a scalar WIMP, we restrict ourselves to 
consideration of unpolarized production of the ^Pi state \E'm (where m = — 1, 0, 1 is the spin index) through 
// fusion at leading order in the velocity expansion. 




where the sum over X is over all SM states in addition to the WIMPonium itself, and the production 
coefficient Djj{^Pi) (which has dimensions of l/[Mass]^) depends on the event kinematics associated with 
the WIMPonium and/or other SM particles involved in the reaction. 
Similar, a fermion WIMP has operators representing its decays, 

1^ <-^(io) + 0{ Si), (16) 

i 

with 

o+cs,) = (e^r)io>(oi(r^c), (i?) 

O-CS,) = (r]=^7y)|0)(0|(r7t^^), (18) 
O°('^o) = (e^^)|0)(0|(7yt^), (19) 
OCS,) = (7?^<?OlO>(0|(?^a^), (20) 

and p-wave state contributions are sub-leading in the velocity expansion. Again, restricting ourselves to 
unpolarized production of a vector WIMPonium state from an // initial state, the relevant operator for 
production is, 



E I*™ + ^) + ^1 (^^'^'?) ■ (21) 



M3 ^ ' ^ 

V X.m 



C. Matrix Elements 



In applying our effective theory to WfMPonium production and decay processes, we encounter matrix 
elements of the non-relativistic WIMP fields sandwiched between WIMP vacua |0) and WIMPonium states 
1^). These matrix elements can be expressed at leading order in the velocity expansion in terms of a 
regularized WIMPonium radial wave function tpr and its derivative tjj'^ at the origin [SJ [T?] , 

{^\OCSom {^\OeS^m^^W?^{a^Mf (22) 
{^\OCP,m A^^(„^M)5 (23) 

where the estimates for the sizes of the matrix elements follow from the velocity expansion [15] , and are 
appropriate for WIMPonium momentum eigenstates normalized to volume. 



IV. IMPLICATIONS AT FUTURE COLLIDERS 



A. Vector WIMPonium 



Vector WIMPonium states may result in theories of complex scalar (as ^Pi states) or Dirac fermion (as 
'^5'i states) dark matter. At colliders, unlike open WIMP production, which leads to signals of missing 
energy, vector WIMPonium signals are more like Z' searches, leading to resonance signals of SM particles. 
WIMPonium decay widths into a final state Y are obtained at leading order in the velocity expansion from 
the imaginary parts of the Wilson coefficients in the four WIMP operators of Equation (10) or (16 1, 



r(* ^ Y) 



2 Im Cy(iPi) 
M4 



r;*V^|*) - 2alM Im CyCPi) 



for a scalar WIMP in a p-wave state, and 

2 Im Cy{^Si) 



2alM Im Cy{^Si) 



(24) 



(25) 



for a Dirac fermion WIMP in an s-wave state. For the case of decays into light SM fermions (m/ <C M), 
the Wilson Coefficient is given by. 



Im C^-^(iPi) = Im Cf,f{^Si) 



Nf 



1 



247r I Ai 

\ ja 



JL 



(26) 



where Nf is the number of colors of /, three for quarks, and one for leptons. These widths are generally 
very small, but nonetheless WIMPonium is still expected to decay promptly on collider time scales. 

Production of an unpolarized ^' at hadron colliders can be driven by qq which at leading order 

involves the same matrix elements as the WIMPonium decays, 

^(99-*) = ^^^^^T.(^\{v^^^)\^rn + X){^„, + X\[^^Vv)\0)-a^DfjCP,) (27) 

X,m 

for a scalar WIMP and 

^(99-^) = ^^M?'^ E(0| \^rn + X){^„, + X\ (CV^) |0) ~ al D^jCS,) (28) 



X,: 



for a fermion WIMP, where at leading order in QCD, 
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where s is the Mandelstam variable corresponding to the partonic center of mass energy squared and we 
have used the fact that to good approximation, ~ 2M. These "partonic" cross sections are convolved 
with the parton distribution functions as usual in order to obtain the hadro-production cross sections. 

With a specific model (or specific choices of the A/s), we can make predictions for production rates and 
decay BRs of vector WIMPonium. For simplicity, we set A/^ = A/^ = A equal for all fcrmions, and switch 
off the interactions with all other Standard Model fields. We can obtain Tevatron limits from Z' searches 
'16] in the language of the c„ and parameters 17] (which are equal under the assumptions that fermions 
all couple universally). For universal fermion couplings and no other decay modes, we have 



,Cd~^a^ (scalar WIMP) c„,Cd-^a^ (fermion WIMP) (30) 

45 45 



where the prefactors assume decays into top quarks are allowed and these can be interpreted as upper limits 
since there are probably fairly large branching fractions into the mediator states whose signatures are model- 
dependent. The CDF limits on c„ and c^i range from a few times 10"^ at resonance masses of 150 GeV to 
about 10^2 at TeV masses. In terms of a^, this requires < 0.2 (0.06) for a scalar (fermion) WIMP whose 
mass is M = 75 GeV and allows essentially any perturbative for WIMPS of mass M = 500 GeV. 

The operators 

{x*'d't^x) (H^Dt'H) and [xi'^x] {H^D^,H) will lead, after electroweak symmetry-breaking, 
to mass mixing between a vector WIMPonium state and the Z boson, with the effective mass^ matrix taking 
the form in the Z-"^ basis, 



(31) 



where n = 5/2 for a scalar WIMP and n = 3/2 for a fermion WIMP. Diagonalizing this matrix (in the limit 
Mz ^ Mxf) leads to Z-^ mixing of the order of a'^M'^/{Ag'^zA?jju). Electroweak precision data generically 
constrains the properties of the Z boson to lie within one part per mil or so of the SM values [TS]. The 
operator [xt'^^'^x] leads to kinetic mixing between $ and the hyper-charge boson, which is much less 
constrained than mass mixing. In the absence of any other WIMP couplings to the Standard Model from 
the UV theory, such mixing may provide the dominant contributions to vector WIMPonium production and 
decay. 

At the LHC, ^' appears like a weakly coupled Z' provided it has appreciable decay into t^t" pairs. The 
cross section will be suppressed by (fermion WIMP) or (scalar WIMP), and thus is expected to be 
much smaller than open WIMP production. However, depending on the underlying theory, the resonance 
in lepton pairs may in fact turn out to be easier to pick out from the background than the missing energy 
signatures of open WIMP production. 



In order to simulate the WIMPonium signal, we reproduce the cross sections at the LHC in equations (27 1 



and (28 1 using MadEvent [IS] to produce ^ ^ ^ and SM background, which we reconstruct with 

the default LHC detector simulation of PCS [20] . An example reconstructed invariant mass distribution of 
for signal and background for a M = 100 GeV Dirac fermion WIMP with universal WIMP couplings 
to fermions Kf — A4 and ^ 0.1 is shown in Figure [2] 

We determine the minimum universal effective coupling of vector WIMPonium to fermions (continuing to 
assume universal Af = M and no other open decay modes) for a 5(t discovery of the WIMPonium against the 
SM background (in an invariant mass cut window of ±10 GeV ~ well above both the intrinsic WIMPonium 
width and the expected experimental leptonic invariant mass resolution - around the assumed mass of the 
WIMPonium) at the LHC with 100 fb~^. We plot our results as a function of the WIMPonium mass in 
Figure [3] where the effective coupling is related to the universal A by, 

geff^a^l'^j^ (scalar WIMP) g^f j a]/'^ (fermion WIMP) (32) 

Once a WIMPonium state is discovered, it can reveal a lot about the underlying theory of dark matter. For 
one thing, the mass of a WIMPonium state decaying into leptons can be measured with great precision, in 
contrast to mass measurements of WIMPs in missing energy events, which are expected to lead to far larger 
uncertainties. Further detailed study of WIMPonium branching ratios and production/ decay distributions 
can reveal many of the details of the pattern of couplings to the SM particles !5T] . The detailed couplings can 
be mapped back into the WIMP-SM coupUngs of equations (|5| or ([6| , and from there into the UV theory. 
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FIG. 2: Sample reconstructed invariant mass of the t~ for the Standard Model (blue) and a vector WIMPonium 
state signal for a mass of 200 GeV, universal A/ = M, and ol^ = 0.2 (red) at the LHC for an integrated luminosity 
of 100 fb"'^. Note that the intrinsic width of the WIMPonium and the experimental resolution are much smaller than 
the bin width. 




200 300 400 500 600 700 800 900 

(GeV) 

FIG. 3: Minimum coupling for a weakly coupled vector WIMPononium state to be visible at 5cr against the SM 
background in the £^£~ channel as a function of the WIMPonium mass, assuming universal coupling to all SM 
fermions, no coupling to SM bosons, and no non-SM decay modes (red curve). Also shown are the CDF limits from 
Z' searches as a function of mass (blue points and interpolating curve). 
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B. Scalar WIMPonium 

Scalar WIMPonium states are very generic. They can occur for any type of scalar or spin 1/2 WIMP. In 
the case of a complex scalar or Dirac fermion WIMP with a scalar mediator particle, there may be several 
scalar states, composed of particle-particle (^'^), particle-antiparticle (^'°), and antiparticle-antiparticle 
(^'~) systems. If there is an exact J7(l)x symmetry, will be stable, which can lead to interesting 
phenomena in cosmology and astrophysics. If U{1)^ is broken, the three states will mix with each other, 
and all three will be unstable. 

Under our assumption that the dark sector does not contribute any new sources of chiral symmetry- 
breaking to the SM fermions, scalar WIMPonium is difficult to produce in the laboratory, because it couples 
very weakly to light fermions. The leading contributions to decay widths are identical for fermion and scalar 
WIMPs, 

r(vI/0 ^Y) = ^^^(^y('So) |(o|^t^|vi,0)|2 ^ Im CrCSo) (33) 

(where for illustration, we write the result for a Dirac fermion state) because the leading matrix elements 
for the ^Sq states are identical in all cases. 

Just as a vector WIMPonium state can mix with the Z boson, scalar WIMPonium can mix with the 
Standard Model Higgs, and have an impact on its couplings. Unless is very large, it is difficult to imagine 
that such mixing would produce an effect large enough to be observable at the LHC |22j, though effects 
could easily be visible if the Higgs is accessible to a TeV scale e+e~ linear collider [23] . 

V. IMPLICATIONS FOR COSMOLOGY AND PARTICLE ASTROPHYSICS 

As alluded to in the introduction, the excesses recently reported by the PAMELA and ATIC experiments 
have inspired exploration of Sommerfeld enhancements in order to motivate the large annihilation cross 
sections required by dark matter interpretations of these signals [U [S] . Such large cross sections are seeming 
at odds with the picture of the WIMP as a thermal relic, and require the effective cross section to be smaller 
at the time of freeze-out than today. This is accomplished by Sommerfeld corrections, which are subdominant 
at the time of freeze-out, when WIMP velocities are on the order of ?; ~ 1/5, but grow at later times when 
WIMP kinetic energies become comparable to the long range potential. This enhancement cuts off for very 
low velocities, once exchanged momenta are of order the mass m of the mediator. There is the potential for 
large effects on the expected rate from indirect DM searches in present day haloes [H [51 [7] . 

In addition to boosting the rate of traditional indirect DM searches, the presence of WIMPonium states 
can also lead to new astrophysical phenomena as well. WIMPs scattering in haloes may form WIMPonium 
by emitting a boson (on or off-shell) with energy on order the binding energy of the WIMPonium. As we 
saw in section |ll] the binding energy for TeV mass WIMPs will typically be of order MeV (up to GeV), 
and can produce e+e~ pairs (perhaps contributing the the INTEGRAL/SPI excess [24j) as well as gamma 
rays. If higher radial modes are available, there could be a whole forest of such lines produced by radiative 
transitions between the various states. Transitions will particularly be favored to direct decay through WIMP 
annihilation if the WIMPonium is produced in a higher angular momentum state [B] . While such lines are 
probably challenging to detect, their detailed spectroscopy encodes the physics of the Yukawa potential 
responsible for binding the WIMPs, which as we saw above determines the number and energy spacings of 
the bound states. 

If the states are stable, some fraction of WIMPs in haloes will probably become bound into them. Since 
these states don't decay through annihilation, they will persist for long times, and can enhance transitions 
between different energy levels. The chemical equilibrium between these stable states and free WIMPs could 
also play a role in the freeze-out process of the WIMP and resultant relic density. Such states could also make 
up some fraction of the WIMPs expected to collect in the Earth and Sun. We leave detailed exploration of 
such phenomena for future work. 

VI. OUTLOOK 

While very simple models of dark matter are appealing, the true dark sector may turn out to be non- 
minimal and involve interesting inter- WIMP dynamics, which may prove important in interpreting signals 
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at future colliders or from astroparticle searches for dark matter. In this article we have explored some of the 
interesting consequences of one such non-minimal dark sector - the possibility that WIMPs may have weakly 
bound states, or WIMPonium. Such bound states are a smoking-gun signature of the dynamics responsible 
for a Sommerfeld enhancement of WIMP annihilation, which has recently received much theoretical atten- 
tion because of the experimental results from PAMELA and ATIC. We have further demonstrated that 
WIMPonium is very interesting in its own right. 

We have constructed effective field theories to describe WIMP interactions with the Standard Model, 
allowing us to be agnostic as the underlying nature of the model of dark matter. We went on to use the 
language of NRQCD (as usually appUed to heavy quarkonia systems) to explore the bounds from LEP and 
the Tevatron, and the prospects for future detection at the LHC. While we have restricted ourselves to 
leading order predictions in the WIMP velocity expansion and order of magnitude estimates for the matrix 
elements of operators, the formalism allows these approximations to be systematically improved. 

The properties of WIMPonium reflect many of the underlying properties of the WIMPs themselves, and 
WIMPonium states can leave their trace at the LHC or by influencing the relic density or signals of indirect 
detection of dark matter. The identity of the dark matter is still a mystery. There is no doubt that unravelling 
that mystery will teach us a lot, and reveal the elegant dynamics of dark matter. 
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